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Abstract 

This report concerns the energy of neutron-matter for densities below 0.15/m -3 and temperatures at 
and below WMeV. Separable NN-interactions are obtained by inverse scattering from the experimental 
phase-shifts with specified momentum cut-offs A. Results of Brueckner-Bloch-DeDomicis as well as finite 
temperature Green's function calculations show independence of cut-off for A > 3. Agreement with the 
low-density virial expansion is found. Results of Hartree-Fock as well as second order calculations show 
considerable A-dependence and the agreement with the virial expansion is lost. The "best" first-order choice 
of A is found to be ~ 2.5/m -1 , which agrees with Vi ow k studies. Reasonable agreement between the second 
order and the Brueckner results is also found for this value of A except at low density. Only 2-body forces 
are used in this study. 

o 
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1 Introduction 

The equation of state (EOS) of neutron matter is of particular interest in astrophysical studies. The mass 
of neutron-stars is related to the EOS at high density ~ 0.5 — 0.7/m~ 3 . Lower densities < 0A5fm~ 3 are 
however of interest in problems related to supernova explosions. The zero temperature Brueckner theory is 
a well tested method at such densities. A finite temperature extension of this theory was given by Bloch 
and De Domicis.[T] It is used here and comparisons are made with the finite temperature Green's function 
method. A main difference from the Brueckner method is that the latter includes the spectral broadening 
self-consistently and the consequences of this broadening has to be investigated. Also, the spectral function 
contains all one-body properties. In that sense it supersedes the Brueckner method. It connects seamlessly 
to the non-equilibrium Green's function quantum transport theory being the stationary solution of this time- 
dependent theory. The spectral functions approach a quasi-particle (Brueckner) limit at low temperatures and 
low density and the energy-integrations that are part of the Green's function method then become numerically 
difficult. The Green's function method is therefore used only at high density and temperature and comparison 
£>. . with the Brueckner-Bloch-DeDomicis method shows satisfactory agreement here. 

There are numerous publications related to the energy of a low density neutron gas, mostly at zero temper- 
ature. The results in one of the earliest [2] already agree satisfactorily with the most recent [31 IU [5] at densities 
low enough for the 1 So-state to be dominant, even though the methods have been different. The finite temper- 
ature equation of state has been studied in detailed calculations by Margueron et al[6] and by Baldo et al[9] 
using Bloch-De Domicis methods. It is the focus of the present work. Our study also relates to recent works 
by Tolos et alj4| using Kohn-Luttinger-Ward methods, and by Horowitz and SchwenkflQ], A main purpose of 
the present study is to illustrate the usefulness of separable interactions that are directly related to scattering 
phase-shifts by inverse scattering to investigate the dependence on cut-offs in momentum-space. The study 
should be considered as exploratory. Any "final" precision calculation of nuclear properties should be based on 
2- 3- (and many-body) forces derived from first principles. 

The potentials are fitted to on-shell scattering data and the ranks of the potentials are the minimum required 
for these fits. If off-shell data are available these can also be fitted by increasing the rank. These data could in 
principle come directly from experiments or from meson-theoretical potentials. In case of the 3 Si — 3 -Di-states 
the deuteron provides such off-shell information. This was used in ref[Tl]. For the problem at hand, neutron- 
matter, this is not relevant. Here the main state is the 1 S , q. For reasons stated below this state is expected to 
be a good candidate for a separable interaction. It was indeed found in ref[TT] that the half-off shell reactance 
matrix elements reproduced the Bonn-B potential data. 

The basic equations used for the inverse scattering problem, the effective interaction in Brueckner theory 
(the G-matrix) and the role of the momentum cut-off is presented in Section 2. A short summary of the Green's 
function method with its effective interaction (the in- medium T- matrix) is given in Section 3. Section 4 presents 
the results of the numerical calculations while a summary of the findings is found in Section 5. 
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2 Separable NN-interaction and Brueckner Theory 



Because the 5-states have poles near E = a reasonable ansatz is to assume that these potentials are 
separable. [12] The inverse scattering method allows the numerical construction of a separable potential that 
reproduces a given set of phase-shifts EXACTLY. This method was used in ref. [11] with application to the 
energy of symmetric nuclear matter. The agreement with Bonn-B results was particularly impressive for the 
1 <SV)-state but also showed surprising agreement for all states except the 3 Pi for which it was found that the 
half-shell reactance matrix differed substantially. The Arndt phases [13] were used then and are also used here 
but at low energy supplemented by the phases obtained from the NN scattering length (a nn = — 18.5/m) and 
effective range (r c — 2.68/m). The effect of low-momentum cut-offs was investigated in several subsequent 
papers jT4l [15l [THJ [17] and the connection with Vi ow k was shown. (See for example Fig. 1 in ref. [17]). Readers 
are referred to these earlier publications for details regarding inverse scattering. Only a short summary and 
some important features of the method are shown below. A rank one separable potential is adequate at low 
density with low relative momenta between the nucleons. But the Sq phases turn repulsive above k c ~ 1.6/to 
necessitating a rank 2 potential as a minimum requirement to fit these phase-shifts and others behaving sim- 
ilarly. Below are only presented the rank one formalism applicable for momentum cut-offs A < k c . For larger 
cut-offs that require a rank 2 potential the equations are modified as in previous work |llj using the Bolsterli 
and MacKenzie method [18]. This is (of course) not a unique method of constructing a potential but as stated 
above it was found to give good agreement with the Bonn and presumably with other realistic potentials. [TT] 
For a rank one separable interaction one has, with the cut-off A indicated explicitly as a parameter 



From inverse scattering one finds: 



where 



< k\V A \\ P >= -v(k; A)v(p; A) (1) 



v 2 (k;A) = ^^sm5(k)\D(k 2 )\ (2) 



D(k ) = exp 



l v f A 



(3) 



where V denotes the principal value and 5(k) is the scattering phase-shift. 

Fig. 1 in ref. [16] shows results of Brueckner calculations with this potential in the lSo-state. There are 
two curves. The uppermost shows the potential energy as a function of A with dispersion-correction, the lower 
without this correction. The density is given by kf — 1.35 /to . The difference between the two curves is 
small. What is also important here is that the energy is constant for A >^ 2.5/m -1 . It is however equally 
important that this is a result of a ladder-summation to all ordersi, Brueckner theory. 

The interaction obtained from the inverse scattering depends on the parameter A as shown already by Fig. 
3 in ref. [17] so to first order in the interaction the energy is not expected to be independent of A. This is 
demonstrated in numerical examples below where the second and all order terms are also considered. 

The A-dependence is further illustrated here by Fig. Q]showing < fc|VA|fc > for two values of A. Although the 
two potentials are notably very different from each other they have in common that they both fit the phase-shifts 
EXACTLY for all k < IfmT 1 . The A = 8/m" 1 potential of course ALSO fits the phases from 2 -> KfmT 1 . 

This means that the diagonal of the reactance matrix defined by 

K{k) = ~ iTiK{k) = ^ ■ tan5 ^/ k ^ 

with 

'*<*>- ^"jf^*"* < 5 » 

is independent of A for k < A. Notice that eq. 2] implies a summation to all orders in V. The in-medium 
interaction defined by the Brueckner G-matrix (or by the in-medium T-matrix defined with Green's functions) 
differs from the reactance matrix K by the Pauli-blocking. Neglecting the dispersion-correction which in our 
case is small (see above), i.e. with no self-energy in the propagator, we have 
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Figure 1: These curves illustrate the dependence of the 1 ^-interaction on A as obtained from the inverse 
sacttering. The maximum value of momentum for each of the two curves is at the two chosen values of A, 2 
and 8/m -1 respectively. 
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with 

'^ p ^l^Ir' ik,) f^ kadk ' (7) 

where P is the center of mass momentum and Q the angle-averaged Pauli-opcrator for pp-ladders, being exact 
with only kinetic energy in the denominator. 

With the Pauli-blocking one has to expect G unlike K to be A-dependent. Calculations show however 
near independence even here as long as A > 2fc/.[Hj This is exemplfied by Fig. [2] showing G-matrix elements 
calculated from the two potentials of Fig. [T] The two curves overlap completely even though the potential are 
quite different. The center-of-mass momentum is here chosen to be zero, which gives the maximum effect of the 
Q-operator. 

Further numerical investigations of the A-dependence as well as of ladder-expansions will be shown below. 
A few observations can be done without reference to detailed calculations. Let the density of nucleons — > 0. 
Then Q — > 1. It is interesting to observe that in this limit 

G —y —An ■ 5{k)/k, 

referred to as the phase-shift approximation pjj] with no dependence on A. (Note the limit gives 8 NOT tanS). 
But it is imporant here that this requires a complete (all order) ladder-summation to be made as in eq. (|6j). 

One might expect the phase-shift approximation to be useful for a low-density neutron-gas. It was however 
already found by Sood and Moszkowski 2] that a proper treatment of the Pauli-blocking at low density, (0.1 < 
kf < 0.5fm~ 1 ) and zero temperature reduces the potential energy/neutron by a factor of ~ 1/2 relative to this 
approximation. The dimcnsionless quantity to consider is in fact kf ■ a s . With a s ~ 20 fm for the neutrons this 
would require kf <C .05, i.e. an extremely low density for the phase-shift approximation to be validQ 

A first order approximation G(k) ~< &|Va|& > would on the other hand be valid only for complete blocking, 
or for a weak potential in which case G ~ V ~ S/k. (With A sufficiently small, i.e. Vi ow t weak enough, one 
may hope to reach this limit.) 

Blocking also increases with density so that a first order approximation may become useful in this case. 
This is consistent with the findings by Moszkowski and Scott. [20] They show the second order contribution 
from their long-ranged part to decrease with increasing density. The long-ranged part in their separation- 
method corresponds to our low-momentum interaction with A ~ 2 — 3/m _1 and it also relates to Vi ow t- |21j 
Numerical agreement between the Moszkowski- Scott method and Vi ow k was shown in ref. [17]. 

It was already emphasised above that eq. © for the G-matrix is a ladder summation to all orders of the 
bare interaction. Below we shall also do a calculation to second order using G -> with 

< k\G^\(P,u;)\p >= v(k)v(p)(l - I G (P,u)) (8) 

Note that with a separable interaction there is no difference in computing effort beteween second and all 
orders. 

Calculations will also be shown below to first order i.e. with 



< k\G (1) \{P,u)\p>=v(k)v{ P ) (9) 

All Brueckner calculations were done neglecting the dispersion correction, i.e. without self-energy insertions 
in particle- and in hole-lines. This approximation is justified by the results of ref. [16] showing this correction 
to be small in the 1 So channel. (See also ref.[3]-) But it is not part of the Green's function calculations. There, 
the calculation of the T-matrix involve integrations over the spectral-functions that by their definition implicitly 
contain the self-energies. 

In the Brueckner case a combination of eqs ^ and ([3]) gives [THl [H] • 

with 

1 [ 2kf 2/r , A ,0(^^)-P,^,/ , kv^k-K) 



lG ^ p) = WrJ " 2 <*» A > V-*« k ' 2dk/ + i^m 

Eqs (fTOlfTTj) have the advantage over eqs (|6][7j) in that the integrand in eq. (fTTj) is zero for k! > 2k f because 
the factor Q(k', P) — V is then equal to zero. Although the two sets of equations are numerically identical the 
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Figure 2: G(k, P = 0) calculated from the two potentials of Fig. [T] There is complete overlap of the two 
curves. The fermi-momentum is here kf — 1.0/m" 1 . 



5 



latter set simplifies the computing greatly. The high momentum component of the interaction is eliminated. 
This new equation also shows that G is independent of scaling t>(fc;A) with a constant factor. Fig. Q] shows 
this to be approximately true and it at least partly explains the result shown by Fig. [2l the independence of 
A. Note however that the above eqs (fpQ|) and (fTTjl are obtained only if assuming the neglect of the dispersion 
correction mentioned above. 

Apart from this modified expression, used here at zero temperature the Brueckner calculation of the energy 
proceeds as described repeatedly in the literature. 

At finite temperatures the Pauli-blocking is modified as shown by Bloch and DeDomicis. Only two-body in- 
teractions were included here. Obviously, 3-body etc and higher order graphs become important with increasing 
density. 



3 Green's Function Method 

Expressions and discussions relating to the Green's function equations for the in-medium T-matrix, self-energies 
and total energy of a fermion-system can be found in numerous publications. [23, 25, 26, 27J Only a summary 
of equations used in this work are shown here. 

Like in the previous section the formalism is only shown for a rank one potential with the modifications for 
the rank 2 done as in ref . [TT] . 

With the separable two-body interaction defined above the in medium T-matrix is given by (the dependence 
on A is suppressed but is implicit) 

where P is the ceneter of mass momentum and 

I(P W )= f fcV W dfc [dcos(8) f— /"^" S (Pi. w '- w ")^2,w")(l-/K-a;")-/(a;")) 



(27r) 3 J w J 2tt J 2tt lo-lo' + ir) 



where 

pi, 2 = — + k 2 ±Pkcos(6). 
The imaginary part of the selfcnergy is given by 

ImZ(p,u) = -^- j^-j ^< P_^|/ m T(|p + k|, W + ^)l E y^>^^')(/(^) + ^ + ^))- 

where b{ui) is the Bose function. The real part i?eS is obtained by the usual dipersion relation and the HF-term 
is added. 

With the self-energies given, the spectral function is calculated and a new T-matrix is obtained. The system 
of equations are iterated until convergence. At each iteration the chemical potential is calculated from the 
nucleon density. The correlated distribution function differs from the fermi- distribution and is given by 

n(p) = ^- f f(u)S{u,p)du (13) 



and the removal energy is 

The total energy is given by 



2tt 



rip) = / wf(u)S(w, P )du (14) 

2im[p) J 

E = ( dp[p 2 /2m - r(p)]n{p) (15) 



(2tt) 3 



4 Numerical Results 

4.1 Energy-Density Results; 1 ^. 

Fig. [3] shows the energy per particle of the neutron-gas as a function of density at temperatures 0,3,4,6 and 
lOMeV. The full lines arc from the Brucckncr-Bloch-DcDomicis calculations. The crosses are from the Green's 
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function calculations. The short full lines are virial results obtained from the work of Horowitz and Schwenk 

The Brueckner calculations were done using standard methods with the cut-off A = 3/m 1 and/or A = 
6/m _1 at the highest density. Increasing A did not change the results while a decrease of A showed noticable 
differences as shown below in section 4.2. It can be remarked that the Brueckner calculations with the separable 
potential and inverse scattering are extremely simple and flexible. Any set of inputted phase-shift data provides 
an output of binding energy in seconds. Dependences on A are easily explored. The momentum-mesh was 
0.01/m- 1 . 

For the Green's function calculations, the set of equations (in medium T-matrix, self-energy S's and spectral 
functions) were iterated until the total energy and chemical potential were stationary. This could require 
anywhere from 5 to 20 iterations depending on the input spectral functions. The in medium T-matrix was 
calculated with the separable potential as defined above and with the momentum cut-off kept constant at 
A = 3/m _1 . The w-integrations were typically from —250 to +250MeV with a mesh of ~ lMeV. 

The Green's function calculations are quite lengthy compared with the Brueckner calculations. They require 
a factor of thousands more computer-time. Although the equations (and the calculations) are rather different 
in the two cases, the results are surprisingly similar. This seems to agree (qualitatively) with zero-temperature 
results of other authors. [25] [27] A bonus is that the spectral functions contain more information on one 
particle properties such as distribution functions and removal energies. Using Thouless criterion[22j allows a 
determination of the BCS critical temperature. This was not investigated in detail but it was determined that 
T c ~ 3MeV at n ~ 0.1/m" 3 . The spectral functions approach a quasi-particle limit at low temperatures and 
low density (See Fig. [9] below) and the energy-integrations then become numerically difficult. This can be 
resolved by resorting to the Extended Quasiparticle Approximation (EQP)[23, 24J or a similar method used by 
Bozek [26] . This was not done here as the Brueckner-Bloch-Domicis method appeared adequate for the present 
and also numerically so much easier. 

4.2 A-dependence and low-order expansion; 1 So 

The dependence on the cut-off in momentum-space was discussed above. Here we show numerical results both 
for the complete Brueckner ladder summation, eq. ([6]), as well as for the second and first order expansions, eqs 
® and ©. 

Fig. 2] shows the energy per particle at T = 6MeV with the full G- matrix, eq. ([BJ. Calculations were done 
for A = 2,3,6 and 8/to" 1 with the Fi gure showing all results converging at the lowest densities. At higher 
densities the A = 2 curve is appreciably higher while for A > 3 one sees convergence. 

Fig. [5] shows second order results, using eq. ((8]). They differ appreciably from the full G- matrix calculations 
as does the first order results shown in Fig. [6] The preferred A showing the best over-all agreement with the 
full G-matrix is A ~ 2.5 for G^ and A ~ 3 for G^ 2 \ Comparing Fig. 2] with Fig. [5] one sees almost perfect 
agreement between the second order and full G-matrix results both for A = 2 and 3/to _1 . Comparison with 
the first order results of Fig. [6] shows also good agreement there for A = 2 but less so for A = 3/to -1 . The force 
weakens if A decreases so that higher order terms then of course become less important. This is in qualitative 
agreement with the rcnormalisation results leading to V\ ow k ■ [I] 

But only the full G-matrix result is capable of showing agreement with the virial expansion, at low density. 

4.3 Energy-Density;including L < 7 

Only the ^o-states were included in the results presented above. At zero temperature that is sufficient for 
densities p < 0.03/m~ 3 .(see e.g.ref.[S|). At higher temperatures even the zero density virial expansion gets 
contribution from the higher partial waves [10] and 3N-forces are increasingly important with density and 
temperature increases. Fig[7] shows results including all partial waves L < 7. Comparison with Fig. [3] shows 
the decrease in energy from higher partial waves increasing with density while it is known that the contribution 
from 3N forces is repulsive and also increasing with density with the two effects partially cancelling each other. 
Comparison with Fig. [3] also confirms the statement made above regarding the contribution of higher partial 
waves to the virial expansion. Fig[7] also shows the A-dependence to increase with temperature. 
The T = result for A = 2.5/m -1 agrees well with the NN-only results shown in refs[4] [28], 
Second order calculations were also done for A = 2.5 with almost perfect agreement with the full Brueckner 
as is to be expected from the results shown in section 4.2. 
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Figure 3: The solid lines show the energy at zero, 0, 3, 4, 6 and lOMeV temperatures respectively from Brueckncr 
ladder summations. The crosses arc from Green's function calculations with selfconsistent spectral functions 
and with in- medium T-matrix ladders to all orders defined by eq. (|12|) . Only ^o-states included here. The 
short lines are low-density virial results from ref . [lOj . 
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Figure 5: Shown are second order results eq.© at T — 6MeV. The uppermost curve is for A = 2fm 
followed by A = 3, 6 and 8/m _1 . 1 S'o-states only. 
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Figure 6: Shown are first order calculations eq.((9]) at T = 6MeV. The uppermost curve is for A = 2fm 
followed by A = 2.5, 3 and 

Note that the A = 8 curve is not shown, ft would be much below the others. 1 S'o-states only. 
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Figure 7: Energy as a function of density including all twO-body states with L < 7. The two top curves are 
for T = lOMeV with A = 2.5/m -1 (upper curve) and A = 5/m _1 respectively The two lower curves are for 
T = with the same values of A. The short line at T = 10 M eV is the virial expansion. 
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4.4 Spectral Functions etc 

The output from Green's function calculations give a more detailed information about the many-body system 
than the corresponding Brueckner quasi-particlc calculations. The additional information is contained in the 
spectral functions. Fig. [3] shows however that the quantity of interest here, the energy per particle as a 
function of density, is practically the same from the Brueckner as from the Green's calculations. The Green's 
function calculations are however computationally much more demanding, in particular at low density and low 
temperature where the quasiparticle picture dominates. The figures [8] and [9] both show spectral functions at a 
temperature T = lOMeV but the first is at a high density while the second is at a low density. One may note 
that these spectral functions for the neutron-gas where the two-body interactions are only in the 1 Sq channel 
are rather different from those also including the 3 Si states with stronger correlations. 

5 Summary 

The energy-density relation for neutron matter at densities ranging from zero to ~ 0.14/m~ 3 and at temper- 
atures ranging from zero to lOMev have been computed. Two-body interactions were obtained from the free 
scattering phase-shifts by inverse scattering with momentum-cut-offs of A = 2 — » 8/m _1 . Three-body interac- 
tions were not included. The many-body calculations were made using Brueckners technique extended to finite 
temperature by Bloch and DeDomicis. Insertions in hole- and particle-lines, i.e. dispersion corrections were 
neglected in this exploratory study as they are small for the states and densities considered here. They should 
of course be included in a " precision" calculation. 

The free scattering interaction is in general A-dependent as illustrated in FigfT] while the in- medium Brueck- 
ner G-matrix in Fig. [2] shows less or no dependence. This was discussed and partly explained above in section 
2. Fig. 2] shows however that when calculating the total energy there is a significant difference between the 
A = 2 result and the A > 3/m -1 results. But there is a nice convergence with increased A. 

Our second order result is noticably A-dependent although less so than the first order with no sign of 
convergence in any case and with considerable difference from the Brueckner results. The value of A with the 
best overall agreement with the Brueckner results is ~ 2.5/m" 1 , which seems to agree with the Vi ow k findings. 
It was also confirmed that for such a small value of A second order and Brueckner agree closely except, as 
remarked below, at low density.. 

We note that there is also a close agreement with Schwenk and Pethick. [5] for T = 0. 

It was found, as did Tolos et al[3] that the higher order terms (in V) are necessary to correctly obtain the 
virial limits at finite temperatures associated with the energy minima in Fig. [3] In fact, it was found that 
the full Brueckner summation was desirable. But only 1 So states are included in this figure. Horowitz and 
Schwenk [TO] find that the virial coefficient 62 increases if higher partial waves are also included, thus further 
enhancing the minima. This increase in 62 becomes larger with increased temperature (Table 1. of ref.[10j) and 
is about 20% at the highest temperature T — WMeV in our Fig. [3] The virial corrections shown in this figure 
are extracted from ref. [10] and include all partial waves. The importance of the higher partial waves in this case 
is in (qualitative) agreement with Fig[7] This is quite different from the situation at T = as shown comparing 
Figs [7] and [3] The So clearly dominates at low temperature and density. 

Comparisons of the Brueckner results with Green's function results were made at T = 6 and lOMeV and 
an agreement within numerical uncertainties was seen. 

A comparison between Brueckner and Green's function calculations for the neutron gas at zero temperature 
was done by Bozek et al[25]. These authors found agreement between the two methods for the densities 
considered here and their results also agree with the present work. At higher densities they found the Green's 
function results much more repulsive than the Brueckner. 

The separable interaction used here may be considered phenomenological. It fits experimental on-shell 
scattering data. Off-shell data are not available. As illustarted in ref. [TTJ it does reproduce the half-off-shcll 
data of the meson theoretically derived Bonn-potentials. As such it may be considered "realistic" for the 1 So- 
state. The separable nature of this interaction is supported by the large scattering length assocated with the 
nearly bound state in the 1 S'o channel. [12] Further comparison with the Bonn-potentials for other states was 
found in ref. [11], an exception being the 3 Pi-state, requiring in this case that the rank of the potential be 
increased. 

The 3N forces were not considered in the present work. They are important and should of course be included 
in the determination of the EOS, but are also less well determined than the 2N forces. A main purpose of the 
present study was to shed some light on the dependence on momentum cut-offs and to show that separable 
potentials derived for inverse scattering are useful for such a study. 
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Figure 8: Spectral functions S(p,ui) at T — lOMeV temperature and a density of 0.11/m 3 . The energy 
(u> — /i =Omega) is in units of fi 2 /m.(/U is chemical potential). The momenta are offset and range from p = 
to the cut-off A = 3/to -1 . The broken curve is at the fermi-surface. Compare this with Fig. [5] which is for a 
lower density. 1 S'o-states only. 
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Figure 9: Similar to Fig. [8] but at a density of 0.017/m 3 . Notice the sharply peaked functions indicating 
the closeness to a quasi-particle limit. This decreses the accuracy of the w-integrations and requires special 
treatment, ^g-states only. 
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